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Superconductivity and magnetism are usually the conflicting (competing) phenomena. We show, 
however, that in nanoscopic objects the electron pairing may promote the magnetic ordering. Such 
situation is possible at low temperatures in the quantum dots placed between the conducting and su¬ 
perconducting reservoirs, where the proximity induced electron pairing cooperates with the correla¬ 
tions enhancing the spin-exchange interactions. The emerging Kondo resonance, which is observable 
in the Andreev conductance, can be significantly enhanced by the coupling to superconducting lead. 

We explain this intriguing tendency within the Anderson impurity model using: the generalized 
Schrieffer-Wolff canonical transformation, the second order perturbative treatment of the Coulomb 
repulsion, and the nonperturbative numerical renormalization group calculations. We also provide 
hints for experimental observability of this phenomenon. 

PACS numbers: 73.23.-b,73.21.La,72.15.Qm,74.45.+c 


Introduction 


Correlated quantum impurity immersed in the Fermi 
sea usually develops the spin-exchange interactions^, that 
cause its total (or partial) screening below some char¬ 
acteristic (Kondo) temperature This effect is 

manifested in the quantum impurity spectrum by the 
Abrikosov-Suhl peak appearing at the Fermi level. It has 
been predicted^ and experimentally confirmed^— that 
in the correlated quantum dot (QD) embedded between 
metallic electrodes, such effect enhances the zero-bias 
tunneling conductance^. This situation changes, how¬ 
ever, if one (or both) external electrode(s) is (are) su¬ 
perconducting because of the proximity induced electron 
pairing^iIS. Depending on the energy level £d, Coulomb 
potential Ud and the coupling Tg to superconducting 
reservoir, the ground state may evolve from the spinful 
configuration | a) (where a =t,f) to the spinless BCS- 
type state u d |0) - v d Such quantum phase tran¬ 

sition (QPT) has a qualitative influence on the spin¬ 
screening mechanism—!. In this work we show that, for 
Ts < Ud , the proximity induced electron pairing can 
strongly amplify the screening effects and give rise to 
a broadening of the Kondo pea k -— » — (see fig. [T| . 

At first glance, such tendency seems to be rather 
counter-intuitive because Tg supports the proximity in¬ 
duced electron pairing that should compete with the 
magnetism. We provide microscopic arguments explain¬ 
ing this intriguing result, based on three independent 
methods. Our study might stimulate and guide future ex¬ 
perimental attempts to verify this theoretical prediction 
in the N-QD-S heterostructures (schematically displayed 
in the left panel of fig.[U), using e.g. self-assembled InAs 
quantum islands^, semiconducting quantum wire s 15 : 16 
or carbon nanotube o 17 ’ 18 . Former measurements of the 
subgap differential conductance have already provided 
evidence for the Andreev/Shiba bound states^ - — and a 



FIG. 1: (Color online) Schematic view of the energy spectrum 
of the N-QD-S junction in the spinfull doublet configuration 
(top panel), where the QD Andreev bound states (driven by 
the coupling Fs to superconducting reservoir) coexist with 
the zero-energy Kondo peak (originating from the Coulomb 
potential Ud and the coupling Tjv to metallic lead). Bottom 
panel shows, that width and height of the Kondo resonance 
strongly depend on Ts. 


tiny (but clear) signature of the zero-bias anomaly driven 
by the Kondo effec t 14 i 16 ' 22 i 23 . Its variation with respect 
to the ratio Ts/LQ has not been investigated carefully 
enough, but this seems to be feasible. 

Similar zero-bias anomalies driven by the super¬ 
conducting proximity effect are nowadays intensively 
explored also in the quantum wires coupled to the 
s-wave superconductors, signaling the Majorana-type 
quasiparticles^— . These exotic quasiparticles originate 
solely from the Andreev/Shiba states in the presence of 
the strong spin-orbit interaction and the Zeeman effect^. 
The present study might hence be useful for distinguish¬ 
ing the zero-bias enhancement due to the Kondo effect 
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from the one driven by the Majorana-type quasiparticles. 


Results 

In what follows we address the proximity induced elec¬ 
tron pairing and study its feedback on the Kondo state, 
focusing on the deep subgap regime. First, we introduce 
the model and discuss its simplified version relevant for 
the deep subgap states. Next, we discuss the issue of 
singlet-doublet quantum phase transition in the limit of 
negligible coupling to the normal lead, r^v —> 0, empha¬ 
sizing its implications for the Kondo-type correlations. 
We then determine the effective spin exchange poten¬ 
tial, generalizing the Schrieffer-Wolff transformation^ for 
the proximized quantum dot, and confront the estimated 
Kondo temperature with the nonperturbative NRG data 
(showing excellent quantitative agreement over the re¬ 
gion Tg < 0.9f7d). We also discuss the results obtained 
from the second-order perturbation theory (SOPT) with 
respect to the Coulomb potential, that provide an inde¬ 
pendent evidence for the Kondo temperature enhance¬ 
ment by increasing Tg (in the doublet state). Finally, we 
discuss the experimentally measurable conductance for 
the subgap regime and give a summary of our results. 


Microscopic model in the subgap regime 

For the description of the N-QD-S junction we use the 
Anderson impurity model? 8 

H = EA+E £ d d\da + Ud hd^hdi 

0 er 

+ d^Ckcrp + ^k/3^ko73 < ^ T ) » W 

k,cr (3 

where /? refers to the normal (j3 = N ) and superconduct¬ 
ing (/3 = S) electrodes, respectively. The operator da 1 
annihilates (creates) an electron with spin a and energy 
e d in the quantum dot, while Vk 0 denotes the tunnel¬ 
ing matrix elements. The repulsive Coulomb potential is 
denoted by Ud and fida = d\d a - Itinerant electrons of 
the metallic reservoir are treated as free fermions, Hn = 
)Ck a ZkNCkaj^CkaN, and the isotropic superconductor is 
described by the BCS model H s = XI k .<7 £kSC ko -sCk<rS - 

Z) k A (4ts fi -kts + c-ktsckts) • Here, denotes the 
annihilation (creation) operator of a spin-er electron with 
momentum k of energy £k /3 in the lead /?, while A 
denotes the superconducting energy gap. It is conve¬ 
nient to introduce the characteristic couplings = 
27rJ] k |I4^| 2 S(co— £k,a), assuming that they are constant 
within the subgap energy regime |w| < A. 

Since we are interested in a relationship between the 
Andreev/Shiba quasiparticles and the Kondo state we 


can simplify the considerations by restricting ourselves 
to an equivalent Hamiltonian^ 

H — Hjsj ^ ) (VkAT d^CkaN T 1 kkT.Y ^ 

k.CT 

+ y t Eddjdg + Ud hd^fidi - (a dd\d[ + h.c.) (2) 

CT 

relevant for the subgap regime in a weak coupling limit 
Ts < A. Effects due to the superconducting electrode 
are here played by the induced on-dot pairing gap A^ = 
r,j/22iiLL. This Hamiltonian © neglects the high-energy 
states existing outside the energy gap window |w| > A 
(see Methods for a discussion) that are irrelevant for the 
present context. 


Subgap quasiparticles of the proximized quantum 
dot 


To understand the influence of electron pairing on the 
Kondo effect, it is useful to recall basic aspects of the 
singlet-doublet quantum phase transition in the ’super¬ 
conducting atomic limit’ Tn —> Exact eigenstates 

of the proximized QD are then represented either by the 
spinfull configurations |cr) with eigenenergy Ed, or the 
spinless (BCS-type) states 

|-) = Ud |0> - Vd in), (3) 

|+) = |0) + u d |n) i (4) 

whose eigenvalues are 

E t = ^Ed + — ^ T \J ^Ed + + a \ , (5) 

with the BCS coefficients 

£ d + Ud /2 

( £ d + Ud/ 2) 2 + A? d 


= l-vj. (6) 



The single particle excitations, between the doublet 
and singlet configurations, give rise to the follow¬ 
ing quasiparticle branches ±Ud /2 ± Ed, where E d = 
y/{£d + U d / 2) 2 + (r s /2) 2 . Two energies ± (U d /2 - E d ) 
can be regarded as the low-energy excitations, whereas 
the other ones (shifted from them by Ud) represent the 
high-energy features. In realistic systems (where U d is 
typically much larger than A) the latter ones usually 
coincide with a continuum formed outside the subgap 

regime^— >21 . 

Diagonal part of the single particle Green’s function 
(for its definition see Methods) is in the subgap regime 
given byii 


G nM = 


a ut 


u+i^-Ed) uj-^-Ed)’ 


(7) 
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FIG. 2: (Color online) The normalized spectral function 
A(uj) = ^YnPcl{ijj) of the half-filled quantum dot obtained 
from the superconducting atomic limit solution (using the 
quasiparticle broadening Fjv = 10 -1 Ts) superposed with 
the Abrikosov-Suhl peak whose width Tk is expressed by 
Eq. 123. The solid/dashed lines correspond to the dou¬ 
blet/singlet ground state configuration and the thick-red 
curve indicates the quantum phase transition at Ts = Ud- 


design the canonical transformation 

k = e § He~ § , ( 8 ) 

which perturbatively eliminates the hybridization term 

v = Xlk.cr (^kJV dJckrrJV + ^kV-L/v^) ■ To simplify the 
notation, we skip the subindex N that unambiguously 
refers to the metallic lead. The terms linear in I4./V can be 
cancelled in the transformed Hamiltonian H by choosing 
the operator S from the following constraint 


H 0 ,S 


= V, 


(9) 


where Hq = H — V. For the Hamiltonian 0 this can be 
satisfied with the anti-hermitian operator S = So — Sq, 
where 


where the partial spectral weight is a = [exp (^) 
+exp (/A)] / [2 exp (^) + exp ( = ^ A ) + exp (-fv)] and 
we set the Boltzmann constant equal to unity, ks = 1. 
The missing amount of the spectral weight 1 — a belongs 
to the high-energy states existing outside the supercon¬ 
ductor gap. At zero temperature, the subgap weight 
changes abruptly from a = 0.5 (when Ed < Ud/2) to 
a = 1 (when Ed > Ud/2). At Ed = Ud/2 the quasipar¬ 
ticle crossing is a signature of the quantum phase tran¬ 
sition from the doublet | a) to the singlet configuration 

| \ 9 . 11.13 

For infinitesimally small coupling one can extend 
the atomic limit solution 0 by imposing the quasiparti¬ 
cle broadening G(oj) —> G(oj + |Fjv). Figure [2] shows the 
normalized spectral function A(oj) = ^T^pd^co), with 
Pd(u>) = — 7 r” 1 ImGn(a;), for the half-filled quantum dot, 
£d = — Ud/2 . On top of these curves we have added the 
Abrikosov-Suhl peak (at to = 0) whose width is given by 
the Kondo temperature, see Eq. m- Upon increas¬ 
ing the ratio Ts/Ud, the Andreev quasiparticle peaks 
move closer and they ultimately merge at the critical 
point Tg = Ud, and simultaneously the Abrikosov-Suhl 
peak gradually broadens all the way up to the QPT. For 
Tg > Ud-, the Andreev peaks drift away from each other 
(see the dashed lines in Fig. 0 and the Kondo feature 
disappears for the reasons discussed in the next subsec¬ 
tion. 


So — EZ EE 7“k 

k «=+,— 

"EE 72,k (cl t d\n^ - cj^nSj.) (10) 

k ck—+,— 

andi 


f d\d c , for a = +, 

\ 1 — dld a for a = —. 


( 11 ) 


The second term of Eq. (flT)ll explicitly differs from the 
standard operator used by Schrieffer and WolfS. From 
the lengthy by straightforward algebra we find that the 
constraint 0 implies the following coefficients y“ k 


7l~k 


7 2 tk 


7i, k 

7 2 'k 


(£k + £d)Vk 

£k _ Ud {£,k + £d) ~ ikd + ^d) ’ 
_AdI4_ 

£k + Ud{£k — Ed) — i^d + A|) 

14 _ A d + 

Ck ~£d Ck-£d 72,k ’ 

A d 4 . 

Ck + £d 7l,k ' 


( 12 ) 

(13) 

(14) 

(15) 


Spin exchange interactions and Kondo temperature 

Adopting the Schrieffer and Wolff approach^ to the 
Hamiltonian 0 of the proximized quantum dot we can 


For A d = 0, the coefficients 7 I? k identically vanish and 
the other ones, given by Eqs. (11211 41) . simplify to the stan¬ 
dard expressions 7 ^ = 14/ (£k -Ud~ Ed) and 7)4 = 
14/ (£k — Ed) of the Schrieffer-Wolff transformation^ 

In the transformed Hamiltonian 
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H — Y. ^kC^Cko- + 2 7l,k^ P (^ka^P^ + l 1 - 0 -) + ( £ d 'Y, Tl.k^k j d^dg 

kcr kpa a \ k / 


A d + Y ^ 2 ,k^k ) (d\d\ + h.c.) + 


k 

w + 


- *e 


7 2 ,k 


v r 


k P 


[(4 


-cl I + cLc£, ) + h.c. 


+ \ v ( c 1 c f 

kn kt‘ p 


+ 2 S (^2,k 72.ky K P 

kp 

+ ^ H ( 7 -tk - 7 r k ) H 

kpa 


C/ d + 2 ^( 7 l -k- 7 l lk)^ 

k 

\U^ - + 

kpcr 

At A it. 


ndt™c4 

k — 7 yk ) 4 (4cVp-<4- CT d<7 + h.C. 


kf p4 ^ pt H 

iA^ +4t4j.44+444t 6 Pt+44 £ h £ p 4+ h - c - 

■ p -<j) 

1 - 


(4(7Cp<T^1(4-0- + C^cr^^-frCp-cr ) + h.C. 


(16) 


we can recognize: the spin exchange term, the interaction 
between QD and itinerant electrons, the pair hopping 
term, and renormalization of the QD energy and the on- 
dot pairing. Since we focus on the screening effects, we 
study in detail only the effective spin-exchange term 

Hexch 1 ( 7kp S d * Skp, (17) 

k.p 


where describes the spin operator of the dot and Sk P 
refers to the spins of itinerant electrons in metallic lead. 
Other contributions are irrelevant for the Kondo physics. 
Formal expression for the effective exchange potential 


7k, P = - 


\ ( 7 i + k - 7 i , k ) V p + ( 7 i + P - 7 i , P ) 4 ( 18 ) 


To estimate the effective Kondo temperature in the 
case of spinfull configuration (for Tg < Ud ), we use the 
formula ! 2 ' 33 TV = ^D exp {— c\> [2p(ep) 7k F k F ]}, where 
p{ef) is the density of states at the Fermi level, D is 
the cut-off energy and the auxiliary function is defined 
as (j)(y ) ~ — 0.5In|j/|. In present case the Kondo 

temperature is expressed by 


t k = V --- 


exp 


Ed ( E d + U d ) + A^ 

r N u d /Tr 


( 21 ) 


with r) being a constant of the order of unity. Influence 
of the on-dot pairing on the Kondo temperature can be 
well illustrated considering the half-filled quantum dot 
case E d = —U d /2■ The spin exchange potential (THH) is 
then given by 


is analogous to the standard Schrieffer-Wolff result^ 
but here we have different coefficients k expressed 
in Eqs. 1121141) . This important aspect generalizes the 
Schrieffer-Wolf potential^ and captures the effects in¬ 
duced by the on-dot pairing. 

In particular, near the Fermi momentum the exchange 
potential (TT51) simplifies to 


7k F ,k F 


Ud |4 F | 2 
Ed (U d + Ed) + A d 


(19) 


It is worthwhile to emphasize that this formula (fTTJl) pre¬ 
cisely reproduces constraint for the quantum phase tran¬ 
sition discussed in the previous section. To prove it, 
we remark that 7k F ,k F changes discontinuously from the 
negative (antiferromagnetic) to the positive (ferromag¬ 
netic) values at £ d ( U d + £ d ) + A^ = 0. Such changeover 
occurs thus at 



which is identical to the QPT constraint = ( U d /2 ) 2 
originally derived in Ref. 11 . 


= -m\v kF \ 2 

■ t/J — (2A d ) 2 


( 22 ) 


and for A^ = 0 it reproduces the standard Schrieffer- 
Wolff result^ 


7, 


N 

kp ,kp 


4|T4 f | 2 

u d 


( 23 ) 


characteristic for the impurity hosted in the metallic 
reservoir. The relative change of 7k F .k F arising from the 
on-dot pairing is 


7k F ,k F _ _ U_d ___4_ ( 24 .') 

44,k F Uj - (2A d ) 2 “ 1 - (Ts/Udf ' 


For the doublet phase (Fg < U d ) the exchange cou¬ 
pling is antiferromagnetic, whereas for the singlet state 
(Tg > U d ) it becomes ferromagnetic. In the latter case, 
however, such ferromagnetic interactions are ineffective 
because the spinless BCS singlet, Ud |0) — v d |!t}> cannot 
be screened. 

The estimated Kondo temperature m increases ver¬ 
sus r 5 , all the way to the critical point at Ts = Ud • 
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FIG. 3: (Color online) (a) The normalized spectral function 
A(io) of the correlated quantum dot obtained by NRG for the 
model Hamiltonian m at half-filling for different values of Eg, 
as indicated. Note the logarithmic energy scale in panel (a), 
(b) The Kondo temperature Tk extracted from the half width 
at half maximum (HWHM) of A(u) from panel (a) (circles) 
and Tk obtained from Eq. m with p = 0.6 (solid line). 
Tk denotes the Kondo temperature in the case of Eg = 0, 
Tk ~ 10 -5 . The parameters are: U d = 0.1 and Tjv = 0.01. 
All parameters are in units of band halfwidth W = 1. 


Such tendency, indicated previously by the NRG data—, 
is solely caused by the quantum phase transition. In a 
vicinity of the QPT the divergent exchange coupling C2l) 
is a typical drawback of the perturbative scheme. Figure 
[3] demonstrates that the formula (12T1) is reliable over the 
broad regime Tg < 0.9 U d . This straightforward conclu¬ 
sion can be practically used by experimentalists. 


Equilibrium transport properties 

We now corroborate the analytical results with accu¬ 
rate numerical renormalization group calculations. 34,35 
In NRG, the logarithmically-discretized conduction band 
is mapped onto a tight binding Hamiltonian with expo¬ 
nentially decaying hopping, £ ra oc A -n / 2 , where A is the 
discretization parameter and n site index. This Hamilto¬ 
nian is diagonalized in an iterative fashion and its eigen- 
spectrum is then used to calculate relevant expectation 
values and correlation functions. In our calculations, 


1 



FIG. 4: (Color online) The normalized spectral function 
A(uT) of correlated quantum dot obtained by NRG calcula¬ 
tions for A > Fs plotted as a function of energy u and Tg. 
The Kondo peak is present in the doublet region, Fg < Ud, 
while in the singlet region, Fg > U d , the Kondo peak no 
longer exists. The parameters are the same as in Fig. [3] 


we assumed A = 2 and kept Nk = 2048 states dur¬ 
ing iteration exploiting Abelian symmetry for the total 
spin zth component.— Moreover, to increase accuracy 
of the spectral data we averaged over TV, = 4 differ¬ 
ent discretizations 1 3 L — We also assumed flat density of 
states, p = 1 / 2 W, with W the band half-width used as 
energy unit W = 1, Ud = 0.1, Tat = 0.01 and zero tem¬ 
perature. In the absence of superconducting correlations, 
T# = 0, this yields the Kondo temperature, T£- ss HP 5 , 
obtained from the half width at half maximum (HWHM) 
of the dot spectral function p d (w) calculated by NRG. 

Figure |3l(a) presents the energy dependence of the nor¬ 
malized spectral function A{oj) of the correlated quantum 
dot at half-filling for the model Hamiltonian ([2}) calcu¬ 
lated for different values of Ts. In the case of Tg = 0, 
A(uj) exhibits Hubbard resonance for u> = ± 1/^/2 and 
the Kondo peak at the Fermi energy, w = 0. It is clearly 
visible that increasing Tg leads to the broadening of the 
Kondo peak. In Fig. [31(b) we compare the relative change 
of the Kondo temperature obtained from the HWHM of 
A(ui) calculated by NRG (circles) and from the approx¬ 
imate formula HID based on the generalized Schrieffer- 
Wolff canonical transformation (solid line). The numer¬ 
ical constant p was estimated to be 77 = 0.6. The agree¬ 
ment is indeed very good and small deviations occur only 
close to Tg = Ud, but then the system is no longer in the 
local moment regime and the Kondo effect disappears. 

The normalized spectral function of the half-filled 
quantum dot in both the doublet, Tg < Ud, and sin¬ 
glet region, Tg > U d , is shown in Fig. [I] In the doublet 
region we clearly observe the zero-energy Kondo peak, 
whose width gradually increases upon increasing rg. Si¬ 
multaneously the Andreev peak (whose width is roughly 
proportional to Tn) moves toward the gap center. In 
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FIG. 5: (Color online) The normalized spectral function A(u) 
at half-filling and T = 0 obtained from the SOPT calculations 
for different values of the Coulomb correlation parameter and 
Fg = Tjv (top panel), and for different ratios of Fs/Tiv with 
Ud = 3Fjv (bottom panel). 


the singlet state, on the other hand, the Kondo peak 
does no longer exist and the Andreev peaks gradually 
depart from each other for increasing T g. The same evo¬ 
lution of the Andreev and the Kondo quasiparticle peaks 
is illustrated in Fig. [2j combining the superconducting 
atomic limit solution with the perturbative estimation of 
the Kondo temperature (l2lT) . 

Broadening of the Kondo peak upon approaching 
the doublet-singlet transition can be independently 
supported by the second-order perturbative treatment 
of the Coulomb interaction term Udh^hdi- The 
first- and second-order contributions have been dis¬ 
cussed in the context of Andree v 39 ! 40 and Josephson 
spectroscopies^—. Here we focus on the Kondo effect, 
studying its evolution near Tg ~ Ud- Diagonal and off- 
diagonal parts of the self-energy can be expressed by^ 


■ r » ■ - [ ImSjfV) gj 


Sll H = 


(25) 

2u(») = ~ + u d (iA) + — J 

2 TT J U> — UJ 1 + l0 + 


with 


(26) 


1 


I m ^U(21)( W ) Hi(w + uj')p^ 2(21 ) (u}') + n 2 (w + u')p 22(21) (u}') 


dui'. 


—oo 
oo 


Hl(2)M = TT / p 11 (+) (w , )p 2 2 (+) (w-W , )-Pl2 (+) (w , )P21 +) ( W - W, ) 


/ 

UJ . 


(27) 

(28) 


In equations (127I28|1 we have introduced 
pfj{u) = -7r _1 ImG^ f (w)/ ± (w), where / ± (w) = 

[1 + exp (±w/T)] _1 denotes the particle/hole Fermi- 
Dirac distribution function and G FF (w) is the 
Green’s function obtained at the Hartree-Fock level 
HhF = H N + YJ (VkNd-tckcrN + ^kJV^ko-iV^) + 


“b UdiP'd—cr)') ^da 


(A d - U d (d±d t)) d\d\ + h.c. 


When calculating the convolutions mm we have 
used the identities E 22 (w) = -[Sn(-w)]* and 

Si 2 (w) = [S 2 i(—w)]*. 


Figure [5] shows the spectral function A(u) obtained 


from the numerical self-consistent solution of Eqs. (1251 - 
1251) . For comparison with the NRG results we focused on 
the half-filled quantum dot (fida) = 0.5. In the weakly 
correlated case Ud < Ts (corresponding to the spinless 
BCS-type ground state) the subgap spectrum is charac¬ 
terized by two Andreev states (shown by the dashed-line 
curves). For Ud ~ Tg, these Andreev states merge, form¬ 
ing a broad structure around the zero energy. In the 
strongly correlated case Ud >Tg (corresponding to the 
spinfull doublet configuration) we observe appearance of 
the Kondo feature (at zero energy) that coexists with the 
Andreev states— We also notice that the width of the 
zero-energy peak (i.e. Tj<) depends on the ratio Ud/Ug 
and such tendency qualitatively agrees with our estima- 
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corresponding to doublet and singlet states. While ap¬ 
proaching the QPT from the doublet side, we observe 
that the zero-bias Kondo peak is gradually enhanced, 
and its width significantly broadens. This tendency is 
caused by the characteristic Kondo temperature, which 
increases with increasing Ts/Ud■ For Ts > Ud , however, 
the Kondo feature is completely absent (in agreement 
with NRG and Schrieffer-Wolff estimations). The mag¬ 
nitude of the subgap Andreev conductance approaches 
then the maximum value 4e 2 //i near the Andreev/Shiba 
states. We notice the quantitative difference between the 
subgap transport properties (shown in Fig. O and the 
electronic spectrum (displayed in Figs. |4]and[5|). Observ¬ 
ability of the Kondo enhancement would be thus possible 
only close to the QPT on the doublet side. 


Discussion 


FIG. 6: (Color online) The bias voltage dependence of the dif¬ 
ferential Andreev conductance Ga obtained from the SOPT 
calculations for the half-filled QD at T = 0. The left panel 
presents Ga for different values of Coulomb correlation and 
Tg = Tat, while the right panel shows the subgap conductance 
for different values of the coupling to superconducting lead T s 
and for Ud = 3Tjv. We notice that the zero-bias feature in¬ 
duced by the Kondo effect is present only for Ts < Ud (in 
the spinfull doublet) and its width gradually broadens with 
increasing Ts/Ud- 

tions based on the Schrieffer-Wolff transformation and 
with the nonperturbative NRG data. 

Differential Andreev conductance 

We now analyze how the observed features reveal in 
the nonlinear response regime. For possible correspon¬ 
dence with the experimentally measurable quantities we 
consider the subgap Andreev current 

UiV) = ^Jdu T A (u) [f + (u,-eV)-f+(tu+eV)] ,(29) 

driven by the applied bias voltage V. The Andreev 
transmittance depends on the off-diagonal (anomalous) 
Green’s function Ta{uj) = |Gi 2 (w)| 2 . We have com¬ 

puted the differential conductance Ga(V) = dI g^' > de¬ 
termining the non-equilibrium transmittance by the tech¬ 
nique described in Methods. 

Figure [5] shows the qualitative changeover of the sub¬ 
gap conductance for representative values of Ud and T#, 


We have studied the influence of the electron pairing 
on the Kondo effect in the strongly correlated quantum 
dot coupled (by Tjv) to the metallic and (by Tg) to su¬ 
perconducting reservoirs by three independent methods. 
The proximity induced on-dot pairing and the Coulomb 
repulsion Ud are responsible for the quantum phase tran¬ 
sition between the (spinless) BCS-like singlet and the 
(spinfull) doublet configurations, depending on the ratio 
of Tg/Ud- Upon approaching this quantum critical point 
from the doublet side, one observes the enhancement of 
the Kondo temperature with increasing We have 

provided the microscopic arguments supporting this be¬ 
havior based on the generalized Schrieffer-Wolff canonical 
transformation. This perturbative treatment of the cou¬ 
pling to metallic lead revealed enhancement of the an¬ 
tiferromagnetic spin-exchange potential, responsible for 
the Kondo resonance. We have compared the estimated 
Kondo temperature with the numerical renormalization 
group calculations, and found excellent agreement over 
the broad regime Tg < 0.9Ud- We have confirmed this 
tendency (for arbitrary r^v) using the second-order per¬ 
turbative treatment of the Coulomb interaction. Our an¬ 
alytical estimation of the Kondo temperature m can be 
quantitatively verified in experimental measurements of 
the differential Andreev conductance. We have shown, 
that the zero-bias enhancement of the subgap conduc¬ 
tance (already reported^ 4 ’ 1 — > 22 > 2 - for some fixed values of 
Ts) would be significantly amplified with increasing ratio 
Ts/Ud , but only on the doublet side. Such behaviour is in 
stark contrast with the zero-bias anomaly caused by the 
Majorana quasiparticles resulting from the topologically 
non-trivial superconductivity. 
























Appendices 


The deep subgap regime |cj| A 


When studying the proximity effect of the Anderson-type Hamil¬ 
tonian CD one has to consider the mixed particle and hole degrees of 
freedom. This can be done, by defining the matrix Green’s function 
G(r,r') = {<^(r);^(r')» in the Nambu spinor representation, 
= (d|,d^), 'fcd = Here we determine its diagonal and 

off-diagonal parts in the equilibrium case (which is also useful for 
description of the transport within the Landauer formalism). The 
Fourier transform of the Green’s function G^(t, t') = — r') 

can be expressed by the Dyson equation 

G_1 (-)=( V (3°) 


The self-energy accounts for the coupling of the quantum 

dot to external reservoirs and for the correlation effects originating 
from the Coulomb repulsion Ud- 

The quantum dot hybridization with the leads can be expressed 
analytically by E^^^w) = where g k p(w) 

are the (Nambu) Greens’ functions of itinerant electrons. In the 
wide-band limit this self-energy is given by the following explicit 
formula 11 ! 40 





x < V iH 


(31) 


Equation iED describes: (i) the proximity induced on-dot pairing 
(via the term proportional to T$) and (ii) the broadening (finite 
life-time) effects. The latter come from the imaginary parts of 
self-energy ED and depend either on both couplings Fp=N,S (f° r 
energies |w| > A) or solely on Fjy (in the subgap regime |cj| < A). 

In the subgap regime |cj| < A the Green’s function of uncorre¬ 
lated quantum dot acquires the BCS-type structure 


G(u) 


Cj + iTjv/2 — £d 

fs /2 


fs /2 \ 

Cj + iFjsj/2 + £d J 


(32) 


with Cj = to + 7 =-: and f s = F g — , A - . The resulting 

z \/A 2 -w 2 \/A 2 -w 2 

spectrum consists of two in-gap peaks, known as the Andreev 11 ! 19 
or Yu-Shiba-Rusinov 20 i 21 quasiparticles. Their splitting is a mea¬ 
sure of the pairing gap A^ induced in the quantum dot. Figure 0 
displays the characteristic energy scales of the uncorrelated quan¬ 
tum dot. 

For infinitesimally weak coupling F jv = 0+ the in-gap states 
have a shape of Dirac delta functions (i.e. represent the long-lived 
quasiparticles). Otherwise, they acquire a finite broadening pro¬ 
portional to Fjy. In the absence of correlations (for Ud = 0) the 
quasiparticle energies Ea,± can be determined by solving the fol¬ 
lowing equation 44 ! 45 


Ea,± + 


(T s /2)E At± 

V A2 ~ E i± 



(r s /2) 2 A2 

A2 - E %±' 


(33) 


In the strong coupling limit, Tg A, we can notice that in¬ 
gap quasiparticles appear close to the superconductor gap edges 
Ea,± — ±A, whereas in the weak coupling limit, Tg < A, they 

approach the asymptotic values, Ea,± — ±-\/£d + (Fs/2) 2 ■ For 
Fn —> 0 , the latter case is known as the ’superconducting atomic 
limit’. The self-energy ED simplifies then to the static value 


E ° (aj) 2 ( ir N ) ’ (34) 

therefore the Hamiltonian CD can be formally modeled by its fully 
equivalent form 0 . describing the proximized quantum dot cou¬ 
pled to the metallic lead. 



-A 0 A 


FIG. 7: (Color online)Spectral function A(to) of the uncorre¬ 
lated QD obtained for £d = 0, Fn = 10 -3 and A = Fs- The 
dashed line shows the reference spectrum in the absence of 
superconducting correlations, A = 0. The in-gap states are 
separated by 2Ad- 


Influence of the correlation effects 

We note that the early studies of the nontrivial relationship 
between the Coulomb repulsion and the proximity induced elec¬ 
tron pairing of the normal metal - quantum dot - supercon¬ 
ductor (N-QD-S) junctions have adopted variety of the theo¬ 
retical methods, such as: slave boson approach 47 ! 48 , equation 
of motion!!, noncrossing approximation!!, iterative perturbation 
technique!!, path integral formulation of the dynamical mean 
field approximation!!, constrained slave boson method!!, numeri¬ 
cal renormalization group!! -13 ! 30 , modified equation of motion!!, 
functional renormalization group!!, expansion around the super¬ 
conducting atomic limit!!, cotunneling treatment of the spin¬ 
ful dot!!, numerical QMC simulations!!, selfconsistent pertur¬ 
bative treatment of the Coulomb repulsion!! and othe»!i!!i!!. 
Amongst them only the numerical renormalization group (NRG) 
calculations!! suggested the Kondo temperature to exponetially 
increase with increasing Fs when approaching the quantum phase 
transition from the doublet side (r^ ~ Ud)- 

The relationship between the proximity induced on-dot pairing 
and the screening effects can be better understood by analyzing 
the superconducting order parameter (d±d^) and the QD magneti¬ 
zation ( S J) = Q. In Fig. [ 8 ]we show their dependence 

on the coupling Ts for several F^/Ud ratios calculated by NRG. 
For finite superconducting energy gap a sign change of the order 
parameter signals the quantum phase transition!!. However, in the 
case of infinite gap considered here, (d^d^) only drops to zero at 
the transition point 11 ! 30 . As clearly seen in the figure, the order 
parameter (d^cLy) increases from 0 to ! around Fs ~ Ud [Fig.fHIaj] 
corresponding to the QPT. Its enhancement is accompanied by the 
suppression of the dot magnetization, which vanishes in the sin¬ 
glet phase, Fs > Ud , (5^) —> 0, see Fig. 0b). Both the increase of 
(d^d- j-) and the decrease of ( S J) indicate the quantum phase transi¬ 
tion at Ts = Ud- Moreover, it can be also seen that the transitions 
present in the above quantities become smeared with increasing the 
coupling to the normal lead r^v- 


Nonlinear charge transport in the subgap regime 

Under non-equilibrium conditions the Andreev transmittance 
Ta(w) = F 2 N |Gi 2(^)| 2 has to be determined using the lesser and 
greater self-energies!! 

E < ’ > (a;) = E < ’ > 0 (a;)+ E < ’ > 2 (o;) 
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FIG. 8: (Color online)(a) The superconducting order param¬ 
eter (dj.dt) and (b) the magnetization {£>%) of the correlated 
quantum dot calculated by NRG for different coupling to nor¬ 
mal lead IV, as indicated. The parameters are the same as 
in Fig. [3] In panel (b) a small external magnetic field B is 
applied to the system, B = 1CP 6 . 


S <( V) 

E>V) 


where 


fr N f+(cj-eV) 0 \ 

V 0 r N f+(u + eV)J 

.( V N f-(w-eV) 0 

V 0 r N f-(ui + eV) 


v-' < C,>2 

^11 





n < ’ > (V + oj)p22 < (u')dbj' 


4’ >2 h 


L A 

2 ?r 


n < ’ > (aj / + (uj')di •' 


sf 2 ’ >2 H 


.A 

2i r 


n < ’ > (V +u)p>J < (u')dt-’ 1 


oo 

n< ’ >M = i J 


and p<’> = G HF ’ r E<’ >0 G HF ’ a , with G HF ’ r( - a '> denoting the 
respective retarded (advanced) Green’s function. The expec¬ 
tation values ^ dlrdcr ^ and have been determined self- 

oo 

consistently from = -^r J (uj) du and = 

— oo 
oo 

J G <2 {u) duj, where the lesser and greater Greens’ functions 

— oo 

obey G < ' > (ui) = G r (ui)E < ’ > (cj)G a (uj), as discussed in Ref—. 
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